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(jj Indirect searches, and in particular rare decays, have proven to be a 

fruitful field to search for New Physics beyond the Standard Model. While 
the down-quark sector (B and K) have been studied in detail, less attention 
was devoted to charm decays due to the smaller expected values and higher 
Q\ theoretical uncertainties of their observables. Recently a renewed interest 

is growing in rare charm searches. In this article we review the current 
\Q experimental status of searches for rare decays in charmed hadrons. While 

00 the Standard Model rates are yet to be reached, current experimental 

limits are already putting constraints on New Physics models. 
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1 Introduction 



The Standard Model (SM) of particle physics is an exceptionally precise theory and 
up to now is able to describe all the available experimental information, with really 
few exceptions. Therefore deviations from this theory belong either to non-explored 
energy regimes or to extremely small effects. While direct searches at current accel- 
erators have not yet led to New Physics (NP) discoveries, hints of new effects might 
hide in rare processes. 

Due to the flavour structure of the SM, flavour changing neutral currents (FCNC) 
are only allowed at loop level, giving the possibility to explore effects where the 
dominant tree level is forbidden. Furthermore they give access, through the virtual 
loops, to higher energy ranges than those accessible to direct searches. FCNC and 
rare decays in general have been extensively studied in the down-type quark sector, 
i.e. in B and K meson decays, and have proven to constrain tightly New Physics 
possibilities. Less attention was devoted in the past to the charm sector, mainly due 
to the fact that SM values of FCNC for charmed hadrons are very small. In fact, the 
absence of a heavy down-type quark which would decouple the short distance effects 
(as the top quark does for the B and K decays) implies that the Glashow-Iliopoulos- 
Maiani (GIM) mechanism pQ is extremely effective for charmed hadrons. On one 
side this leads to the advantage that NP effects can be orders of magnitude larger 
than the SM ones; however it also implies that SM predictions can be dominated 
by long distance contributions due to the propagation of light quarks intermediate 
states. Unfortunately long distance contributions, being non perturbative, carry large 
theoretical uncertainties. 

Many New Physics scenarios can enhance these processes, like Supersymmetry 
(with or without -R-parity violation) [2], Littlest Higgs model [31 H] and Lepto- 
quarks [5]. Furthermore some models lead to enhancements only in the up sector, 
and therefore can only be probed in charm decays. 

As the experimental results improve a renewed interest in the charm sector has 
grown, also due to the recent measurements of CP violation in D decays by the 
LHCb [6] and CDF experiments [7], stimulating the theoretical work. 

We will review in the following the current status of experimental searches for 
rare charm decays. In particular in §[2] we will treat the FCNC c — > u£ + £~ ^ transi- 
tions for charged (£ |2.1[ ) and neutral charm hadrons (§ |2.2[ ). Radiative processes will 
not be treated in the following as the high theoretical uncertainties make them less 
interesting [2], though recently it has been proposed that their study could shed light 
on charm CP violation [8] . In sections [3] and [4] we will describe the di-photon and 
di-lepton D° decays. Finally a brief overview of the current searches for SM forbidden 
decays will be given in section [5] 

*Here and in the rest of this article we always imply charge conjugation unless otherwise stated. 
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As far as the future prospects are concerned, we will not discuss them in this 
article as they are subject of another report at this same conference. 

2 c — > ui + i~ transitions 

We start from the c — > u£ + £~ transitions as they are the most probable to be observed 
at currently running experiments. Short distance contributions for transitions of the 
type c — > u£ + £~ are heavily suppressed by the GIM mechanism. In particular, the 
inclusive branching fractions calculated for the mesonic modes to electrons are [2J: 

B m ^4 e+6 - ^ 2 • 1(T 8 B D o^ x o e+e - ~ 8 • 1(T 9 , 

which would be out of present experimental reach. However, as already said, short 
distance contributions are shaded by the long distance ones. Long distance (LD) 
processes can be represented in the form D — > XV — > X££ where V — <fi, p, oj is 
the intermediate resonance. Typical exclusive branching fractions of the resonant 
processes are at the level of (9(10 -6 ) j2j|9]. As an example the branching fraction for 
D + — > ii + e + e~ is B ~ 2 x 1CT 6 dominated by the 7r + intermediate state. 

However long distance processes, being non-perturbative in nature, carry also large 
theoretical uncertainties so that they make more difficult the comparison between SM 
and New Physics predictions and between these and the experimental results. 

Different New Physics scenarios can enhance c — > u£ + £~ transitions. Already the 
Minimal Supersymmetric Standard Model (MSSM) with R-parity conservation could 
give contributions to c — > u£ + £~ processes. These would come from one-loop diagrams 
with gluino and squarks in the loop [2J [III] and would increase the decay rate at low 
di-lepton invariant mass. However when including bounds from D° — D oscillations 
these enhancements are reduced. Furthermore when considering the hadronic decay 
D — > X£ + £~ the rate is decreased by a factor when A is a pseudoscalar meson 
and hidden by long distance contributions when A is a vector meson [11] . 

On the other hand the MSSM with R-parity violation (ftp) is not so constrained 
and can give large contributions via a tree-level exchange of squarks. Enhancement of 
a factor ~ 3 are possible, for example, in the D + — > decay [UJ. In particular, 

as shown in Fig. [lja), at low and high di-muon invariant mass (or q 2 = m^), i-e. far 
from the resonances, there is a possible high sensitivity to the MSSM^? P . 

Similar increases in the branching fraction can appear also in the D — > pfx + fi~ 
channel; in this channel in addition to the branching fraction also angular observables, 
such as the forward-backward asymmetry for the leptons are thought to be sensitive 
to NP effect, even if only reachable at high yields j2j ITT] . 

Other possible contributions to c 4 u£ + £~ can come from new physics lepto- 
quarks. In particular it has been shown that the D + — > n + p + yT decay is sensitive 
also to couplings which cannot be limited by constraints on the branching fraction 
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Figure 1: Differential rate as a function of invariant mass squared of the two muons 
in the D + — y decay for the SM long distance (LD) prediction, and compared 

with (a) MSSM^? P [llj and (b) Leptoquarks new physics model. The Leptoquark 
contribution is saturated to the experimental limit on the total branching fraction [12] . 



of D° — y [5]. In Fig. [T|(b) the differential branching fraction of D + — y 7i + fi + jj,~ 

is shown in the hypothesis of saturating its experimental limit with Leptoquarks 
coupling, showing that large room for a distinction between this and the LD SM 
contribution is still present at large and low q 2 . 

In the next two sections we will review the experimental searches for c — y u£ + £~ 
transitions divided for simplicity into charged and neutral hadron decays. 



2.1 Searches for Xf -> h^Pf 

The current limits on rare decays of charged charm hadrons are nicely summarized 
by the Heavy Flavour Averaging Group (HFAG) [T3] plots shown in Fig. [2j 

None of the listed channels have been observed up to now and the limits on the 
decays branching fractions are at the level of 10~ 6 or higher. With few exceptions 
the best limit for all of these decay modes come from a comprehensive analysis of the 
BaBar experiment, described in Ref. [14], which we review in the following. 

2.1.1 Search for Xf -»■ /i ± £+£~ at BaBar 

A search for rare decays in the form Xf — y h ± i + t~ was conducted by the BaBar 
collaboration [T4] exploiting data collected at the PEP-II at SLAC. A dataset of 
347 fb _1 of e + e~ collisions at the T(4S') and of 37 fb _1 just below it was used. We 
will not describe the BaBar detector and the reader should refer to Ref. [T5] . 
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Figure 2: Summary of the upper limits on charged charm hadron rare decays compiled 
by the HFAG group. [13] 



Pairs of tracks identified as leptons and a track identified as a charged hadron 
(pion, kaon or proton) were used to build charm hadron candidates. The candi- 
date momentum in the e + e~ center of mass (p*) was required to be larger than 2.5 
GeV/c 2 to remove a large combinatorial background in those regimes. QED back- 
ground was rejected requiring a multiplicity larger than 5. The residual combinatorial 
background, mainly from charm and beauty semi-leptonic decays, was reduced with 
cuts on the \ 2 °f the vertex fit and on the distance of closest approach between the 
leptons. Finally the <fi region was rejected from the invariant mass of the di-lepton in 
order to exclude the resonant (and long distance) part of the decay. After this first 
selection a likelihood ratio was used to further discriminate the signal from combina- 
torial background. The likelihood was built from information on the p* of charmed 
candidate, total energy in the event and flight length significance. The PDFs for the 
signal were built using Monte Carlo (MC) simulations while for the background data 
events in the invariant mass sidebands were used. The cut value on the likelihood was 
optimized independently for each decay mode in order to provide the lowest expected 
upper limit. 

Extended, unbinned maximum-likelihood fits to the invariant mass distributions 
were used in order to extract the number of signal events. The signal pdf was 
parametrized as a Crystal Ball function while the combinatorial background as a first 
order polynomial. For some channels an additional component to take into account 
background from mis-identified non leptonic charm decays was included. In Fig. [3] 
the invariant mass distribution of some of the studied channels is shown together with 
the fit. 

The signal yields were normalized to known hadronic charm decays in order to 
compute the branching fractions: mesons channels were normalized to the Df g s — > 
<j)(—> KK)ir + channel while A+ — > pK~n + decay was used for the A+ decays. 

No significant signal was observed in any of the channels. A 2.6a fluctuation 
was seen in the A+ — > pfi + fi~~ decay which however has a 25% probability to occurr 
when looking at 35 different channels. Therefore upper limits were put on the decay 
branching fractions using a Bayesian approach with flat prior for the event yield in 
the physical region. The obtained upper limits on the branching fractions at 90% CL 
are between 1 • 10~ 6 and 44 • 10~ 6 depending on the decay mode and are the best 
limits up to now for most of the channels. 

2.1.2 D+ -> at D0 

One of the channels in which the best limit does not come from BaBar is the 
D + 7i + fi + {i~ decay were a strongest limit was put by the D0 collaboration [12] . 
This search was conducted with 1.3 ftT 1 of pp collisions at y/s = 1.96 TeV collected 
at Tevatron with the D0 detector. 

Tracks identified as muons were used as starting point after a cut on the transverse 



5 




M(peV)'[GeV/c a ] " MfpuVHGeV/c 2 ] 



Figure 3: Invariant mass distributions for — > h ± £ + £ candidates in the analysis 
performed by BaBar [14J. 
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momentum (2 GeV/c) and on the momentum (3 GeV/c). Two oppositely charged 
muons were combined requiring to form a well reconstructed vertex and to have an 
invariant mass lower than 2 GeV/c 2 . The invariant mass region of the di-muon 
system was excluded from the search of the non-resonant D + — > system 
and used to study the D + — > (f>7i + decay, with (f> — > The di-muon system 

was combined with an additional track passing requirements on its impact parameter 
significance and on its transverse momentum, and the three body invariant mass was 
required to be in the 1.4 - 2.4 GeV/c 2 range. Further selection criteria were applied 
based on the significance of the transverse flight distance of the D + candidate and 
on the collinearity angle, i.e. the angle between the D + momentum and the line 
between primary and secondary vertices. After this selection the n + fi + fi~ invariant 
mass for events in the <p di-muon mass region is shown in Fig. |4](a), where the two 
signals of D + and — > </)tt + can be seen. The D + — > <p7i + — > branching 
fraction was measured to be (1.8 ± 0.5(stat) ± 0.6(syst)) x 1CT 6 . Before considering 
the non-resonant part of the di-muon invariant mass, tighter cuts were applied based 
on the already mentioned variables and on the pion transverse momentum and the 
isolation of the D + \ X = po/ Y2Pcone where the sum of the momenta is extended to 
a cone around the D + candidate. The criteria were optimized in order to obtain the 
best upper limit on the branching fraction. The invariant mass distribution of the 
remaining events is shown Fig. |4|b). No signal above the combinatorial background 
was observed and the following upper limit on the branching fraction was estimated: 
B(D + n+fi+fi-) < 3.9 x 1(T 6 at 90% CL. 

While this limit is still more than 2 orders of magnitude from the Standard 
Model expected rate it constraints already the parameter space of R-parity violating 
SUSY [2J; on the other hand the little Higgs models are an order of magnitude below 
this level and therefore are still not constrained by this measurement [I]. 

2.2 Searches for X c ° ->> h°£+£- 

Results on rare decays from D° decays are all at least 10- years old, with the exception 
of di-lepton and di-photon decays which will be treated separately. A summary of the 
upper limits on the various decay modes is shown in Fig.[5]from the HFAG group [13"] . 
With few exceptions, all the limits on D° — > Xli decays come from two analyses by 
the E791 and CLEOII collaborations which will be discussed in the following. 

2.2.1 Search for D° XU decays at E791 

The E791 collaboration reported in 2001 a search for decays in the forms D° — > Vt r l~ , 
with V = p°, K*°, (p and D° hM + £~, with h = 7r,K p] based on 2 x 10 10 events 
of a 400 GeV/c 7r~ beam on fixed target. Events with secondary vertex well sepa- 
rated from the primary one and from the detector material were selected. A blind 
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Figure 4: Invariant mass distribution of Dfs — > 7i + fi + fi candidates for events with 
the di-muon (a) in the mass region and (b) outside of it [T2] . 
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Figure 5: Summary of upper limits on branching fraction of FCNC D° decays [13J. 
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analysis was optimized exploiting MC simulations of the signal channels and using 
background from the invariant mass sidebands. All the channels were normalized 
to topologically similar hadronic decays to which the very same selection, with the 
exception of particle identification, was applied. MC simulations were used in order 
to calculate the signal and normalization selection efficiencies, while the muon and 
electron identification efficiencies were measured on data. After unblinding the signal 
invariant mass regions no significant signal was found in any of the channels, as can be 
seen from the invariant mass distributions shown in Fig. |6^a). After subtracting the 
combinatorial background, estimated from the sidebands, and the mis-identification 
background, upper limits on the number of signal events were calculated with the 
Feldman-Cousins method and translated into branching fractions with the normali- 
sation channels. 

The upper limits on the branching fractions vary between 3.0 and 55.3 xlO -5 at 
90 % CL depending on the channel, most of which are the best results to date. 

2.2.2 Search for D° ->■ XU decays at CLEOII 

A similar study was conducted by the CLEO collaboration with the CLEOII detector 
at CESR [T7j. A dataset of 3.87 ftT 1 of e + e _ collisions at the T(45) and just below it 
was used to search for events in the form D° — > Xll where X = 7r°, K s , rj, p°, u, K*°, cf). 
In order to build the signal candidates, charged tracks were required to come from 
the primary interaction, with the exception of the ones for the K s . Particle identifi- 
cation of pions and kaons was based on dE/dx and time-of-flight information, while 
calorimeter and muon chambers information was additionally exploited for the lepton 
candidates. Electrons coming from photon and 7r° decays were rejected. The various 
employed mesons were reconstructed in the following final states: K s — > 7r + 7r~, 
p° — > 7r + 7r~, to — > 7T + 7r~7r°, K — > K~tt + , (j) — > K + K~ , and 7r° and r] into 77. 
The so constructed D° candidate was required to come from a D* + — > D°tt + decay, 
in order to reduce the combinatorial background, imposing a cut on the D* + — D° 
mass difference (AM) at 2.0 MeV/c 2 from the nominal value. Cuts on the scaled 
momentum of the D* (po*/ \j E^ eam — M|,* ) and on the angle between the daughter 
tracks and the D° momentum were also applied. 

After this selection the invariant mass distributions of the studied decay modes 
are shown in Fig. |6](b). No signals were observed over the small expected combinato- 
rial background and a 90% poissonian upper limit was put on the number of signal 
events. This limit was translated in a branching fraction by normalizing to the total 
number of produced D°, measured in data in the K~n + final state, and correcting 
for the efficiency of each channel, obtained from simulations. The upper limits on 
the branching fractions of the various modes are in the range between few in 10 -5 to 
10 -4 depending on the channels and are the best limits to date for most of them. 
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Figure 6: Invariant mass distributions for (a) the search of signals in the form 
D° — » Xll performed by the E791 collaboration [TB] and (b) candidates in the form 
D° — > Xll in the search performed by the CLEO collaboration [IT] . 
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3 IT -> 77 



The Standard Model prediction for the short-distance contribution to the D° — > 77 
decay leads to a branching fraction of 3 x 1CT 11 , which however is completely shaded 
by the long distance contributions [10]. Calculations done for the Vector Meson 
Dominance mechanism give a long distance contribution to the branching fraction of 
3. 5^ 6 • 10~ 8 , and other calculations which exploit the Heavy Quark Chiral Pertur- 
bation Theory lead to similar results pIU [18] . 

In various new physics scenarios this process could be enhanced and it has been 
pointed out that within the MSSM scenario gluino exchange could enhance the SM 
rate of more than 2 orders of magnitude [19J. 

3.1 Search for D° — > 77 at BaBar 

The present best limit on the D° — > 77 branching fraction is due to an analysis 
performed at BaBar and reported in Ref. [201 • This analysis is based on 470.5 fb" 1 
of e + e~ collisions at -y/s = 10.58 GeV and 10.54 GeV. The search for D° — > 77 was 
done in parallel with a measurement of the D° — > 7r°7r° branching fraction which is 
also the main background for the di-photon final state. As normalisation channel was 
chosen the D° — > K^n decay for its well known branching fraction. 

MC simulations were used to optimize the search for D° — > 77 decays. D° — > 77 
(D° — > 7T°7r ) candidates were formed with pairs of photons (pions) required to have 
an invariant mass between 1.7 and 2.1 GeV/c 2 (1.65 and 2.05 GeV/c 2 ). Photons 
were required to have a CM energy between 0.74 and 4 GeV. The D° candidates 
were required to come from a D* + — > D°tt + decay and therefore kinematically fitted, 
together with a tt + , to a common vertex constrained to be within the beam spot. 

Banckground from B decays was rejected with a 2.85 (2.4) GeV/c momentum 
threshold on the D* candidate for D° — > 77 (D° — > 7r°7r°), while the large QED 
background was removed with a cut on the tracks and neutrals multiplicity of the 
events. Finally the D° — > 7r°7r° background to D° — > 77 was rejected with an ad-hoc 
7r°-veto which required the photons not to be consistent with coming from a 7r° decay 
after being combined with other photons of the same event. 

The invariant mass of D° — > 7r°7r° candidates is shown in Fig. [7](a); the superim- 
posed fitted PDF is composed of a 3 rd order Chebychev polynomial for the background 
and of a signal PDF which is the sum of a Crystal Ball, a gaussian and a bifurcated 
gaussian function. A signal yield of about 26k events was estimated which led to 
a branching fraction measurement of: B(D° — > 7r°7r°) = (8.4 ± 0.1 ± 0.3) x 10 -4 . 
On the other hand a negative signal yield, consistent with zero, was estimated from 
the fit to the D° — > 77 invariant mass distribution, shown in Fig. [7^b). The mea- 
sured upper limit on the branching fraction of this decay was estimated to be: 
B(D° -)■ 77) < 2.2 x 10~ 6 at 90% confidence level. While this result is not yet 
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Figure 7: Invariant mass distribution for (a) D° — > 7r°7r° and (b) D° — > 77 candidates 
in the analysis performed by the BaBar collaboration [20] . 



approaching the SM rate, it limits the branching fraction to, at most, ~ 70 times the 
SM levels, tightly constraining New Physics models. 

4 D° -> £ + £~ 

While the two-body leptonic signature of the D° — > £ + £~ decays is more appealing 
from the experimental point of view, theoretically they are foreseen to be even more 
rare than c —> u£ + £~ processes, due to the helicity suppression of the phase space. The 
Standard Model short distance contribution for D° — > n~ is at the level of 10~ 18 
and even less for D° — > e + e~. The long distance contributions are relevant in the 
SM branching fraction and in particular the dominant one is due to the two-photon 
intermediate state and the following expression can be derived for the di-muon final 
state [2]: 

S(D°^/iV)-27x lO-^GD ^ 77) (1) 

which leads to an estimate of a minimum theoretical contribution to the branch- 
ing fraction of B(D° — > > 10~ 13 . However equation ([T]) could also be ex- 
ploited to estimate a rough upper limit to this contribution: together with the 
experimental limit on D° — > 77 set by BaBar [20] described earlier, this gives 
B(D° -> /i + /0 77 <6x 10- 11 . 

Many NP models can instead enhance the D° — > £ + £~ branching fraction by 
various orders of magnitude. Moreover it has been noted that the D° — > decay 
can be related directly to the couplings involved in the D° — D mixing [21] so that 
limits on one of the two processes can be used to constrain evaluations of the other. 



12 



As an example one of the largest contributions is from ft p SUSY and gives an estimate 
of [21]: 

„» t , /300 GevV „ /300 GevV 

* 4.8 x 10- 7 x D l-^—j < 4 ' 8 x 10 ( J ( 2 ) 

which is dependent on the mass of the down type quarks super partners. 

Finally, also contributions from Leptoquarks were found to be relevant, and pre- 
dicted branching fractions at the level of 5 x 1CT 7 [22J and are being already strongly 
constrained by experimental limits. 

In the following we will present a preliminary result of the LHCb experiment 
on the D° — > ji + decay, which is the strongest limit up to date, and the results 
obtained by the Belle collaboration on all the D° — > £ + £~ decays. 

4.1 Search for D° fi+fT at LHCb 

The search for the D° — > fi + decay at LHCb is based on a dataset of 0.9 fb _1 of 
pp collisions in LHC at y/s = 7 TeV [23] . An additional sample of about 79 pb _1 was 
also used in order to optimize the selection but not used in the final analysis. 

Pairs of opposite side muons were selected to form D° candidates, which were 
required to come from a D* + — > D°tt + decay. After a first selection based on stan- 
dard geometrical and kinematic variables, a Boosted Decision Tree (BDT) was built 
in order to reject the combinatorial background, combining information from the fol- 
lowing variables: x 2 °f the impact parameter of the D° candidate and pointing angle 
of the D° with respect to the primary vertex, minimum transverse momentum of the 
two muons, minimum x 2 of impact parameter of the two muons, angle of the positive 
muon in the D° candidate rest frame with respect to the D° flight direction. In Fig. [8] 
is shown the distribution of the BDT output for signal MC events and for events in 
the di-muon invariant mass sidebands. The final cut value on the BDT was chosen in 
order to have the best expected upper limit on the D° — > fi + fi~ branching fraction. A 
large background to this decay is the D° — > 7r + 7r~ with both pions mis-identified into 
muons. This background was estimated by measuring the mis-ID probability directly 
on data D° — > K~7t + decays. The double mis-ID probability was estimated to be 
p(D° ->- 7T+7T- ->- nn) = (27.3 ± 3.4 ± 2.0) • 10~ 6 . The D* + ->- D°(-^ ii + it-)tt + decay 
was also exploited as normalisation channel, after being selected with the same selec- 
tion as the signal, with the exception of the muon identification. The D° — > 7r + 7r~ 
yield was estimated with a 2-dimensional fit to the distribution in the AM versus in- 
variant mass plane. A projection of this distribution together with the fit is shown in 
Fig. ^b). The signal and normalisation channel efficiencies were estimated from MC 
simulations and corrected for data-MC discrepancies. The estimate of the D° — > fi + n~ 
signal yield was also done in the 2D AM-M D o plane. The data distribution was fitted 
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Figure 8: In (a) the distribution of the output of the BDT for signal MC simulations 
(continuous line) and combinatorial background from data (crosses) are shown as 
obtained in the LHCb experiment. In (b) is shown the invariant mass distribution 
for D° — > 7r + 7r~ candidates in data |23|. 



with a PDF built upon the following components: the signal was parametrised as a 
gaussian (in M) times a double gaussian (AM), the D° — > 7t + tt~ mis-ID as Crystal 
Ball (M) times a double gaussian (AM), the D° — > K~ir + tail mis-ID as a gaussian in 
both M and AM and the combinatorial background as an exponential (in M) times 
the following function: /(AM) = {^f (1 - e^^ 1 ) + b- (^ - 1). The projection 
of this fit along the D° candidate invariant mass is shown in Fig. [9^a). As no signifi- 
cant signal was observed, the upper limit to the D° — > branching fraction was 
estimated with the CLs method to be: B(D° ->- < 1.3(1.1) x 1(T 8 at 95 (90)% 
CL. In Fig. [9^b) is shown the CLs value as a function of the D° — > branch- 
ing fraction and the expected values and bands. This limit is still several orders of 
magnitude far from the SM predictions but constraints various NP models. 



4.2 Search for D° -> £ + £~ at Belle 

Another search for di-leptonic D° decays was done by the Belle collaboration and re- 
ported in Ref. [23]. This analysis was based on 660 fb~ of e + e~ collisions at the T(45") 
and just below it, collected with the Belle detector at KEKB. D° candidates were built 
from two opposite side leptons in the D° — > , D° — > e + e~~ and D° — > e^fi 11 final 
states. Only D mesons from e + e~ — > cc were used, excluding those from B decays 
because of the higher combinatorial background; these were removed by requiring a 
maximum allowed missing energy which is high in events with B semileptonic de- 
cays due to missing neutrinos. The D° — > 7t + tt~ final state was also reconstructed 
as normalisation and control channel. Standard particle identification criteria, based 
on dE/dx, TOF, Cherenkov light and muon detectors information, were exploited in 
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Figure 9: (a) Invariant mass distribution for D° — > ji + jj," candidates, projected from 
the Mjj, AM plane, and the fit to the distribution: the black line is the total fit func- 
tion composed of the combinatorial background (light grey dashed line), D° — > tt + tt~ 
mis-identified background (dark dashed line), D° — > K~tt + mis-identified background 
(dash-dotted line) and the signal D° — > ^ (light grey continuous line). In (b) the 
CLs value as a function of the D° — > n + fi~ branching fraction as obtained with the 
CLs method is shown; the red line shows the upper limit to the branching fraction at 
95 % CL [23]. 
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order to select pions, electrons and muons. The D° candidate was required to come 
from a D* + — > D°tt + decay from the vertex of the e + e~ interaction point. Further 
rejection of combinatorial background and of B-decays was achieved by requiring a 
D* + candidate momentum in excess of 2.5 GeV/c. Candidates were selected to be 
in a region of the D° candidate invariant mass between 1.81 and 1.91 GeV/c 2 and 
with q < 20 MeV, where q = (Md*+ — M^o — m n )c 2 . Further selection criteria, based 
on the missing energy, lepton identification and signal region in the mass versus q 
plane, were optimized differently for each final state with a blind analysis. Residual 
combinatorial background was estimated from the sideband region of the M — q plane, 
while the mis-identified D° — > 7r + 7r~ background was evaluated from the decay itself 
reconstructed in data, re-weighted for momentum dependent mis-identification prob- 
abilities. The invariant mass distributions after the final selection criteria are shown 



in Fig. 10 No events in excess of the estimated background are observed in any of 
the final states. A binned maximum likelihood fit to the invariant mass distribution 
was used to estimate the yield of the normalisation channel D° — > 7r + 7r~. Efficien- 
cies were estimated from tuned MC samples. Finally the following upper limits were 
derived on the signal channels branching fractions: B(D° — > < 1.4 x 10 -7 , 

B(D° e + e~) < 7.9 x 10~ 7 , B(D° -> e ± ^) < 2.6 x 10~ 7 at 90% CL. 

The ones on D° — > e + er and D° — > e ± /i =F are the best limits on these channels to 
date and limit further the parameter space for NP models. 



5 Forbidden decays 

While rare decays are often dominated by the theoretical uncertainties, an unambigu- 
ous sign of New Physics would be to observe one of the decays which in the Standard 
Model are forbidden because of violation of one or more conservation laws. In partic- 
ular here we consider decays which do not conserve the Lepton Family Number, the 
total Lepton Number, the Baryon Number or more than one of these. 

Experimentally none of these decays has ever been observed, and therefore limits 
on the branching fractions are set. A summary of the current situation for the neutral 



sector, i.e. for D° decays, can be seen in Fig 11 where is shown a collection of the 
various limits on branching fractions compiled by the HFAG group [13]. These upper 
limits have been set, with few exceptions, by the already mentioned experiments: 
E791 HB1 and CLEOII 



Analyses of D — > X££ LFV decays were performed in the very same way as 



described in £2.2.1 and ^2.2.2 for the E791 and CLEOII experiments respectively. 



The invariant mass distributions for some of the analysed final states are shown in 



Fig. 13 None of the searched decays was observed and the limits on the branching 
fractions are between 10~ 5 and few in 10~ 4 . 

One of the exceptions in this sector is the already mentioned limit on D° — > e fi^ 
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Figure 10: Invariant mass distributions for D° — > £ + £ candidates in the analysis 
performed by the Belle Collaboration [24] . 
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Figure 11: Summary of branching fraction upper limits on D° forbidden decays, 
compiled by the HFAG group [13]. 



which has been set by the Belle collaboration [21] and described in £4.2 

The BaBar collaboration has instead the almost absolute monopoly on the limits 
on the charged sector as can be seen by the HFAG summary plots in Fig. [2j Most 
of the limits on LFV and BV decays for D + ,Df and A+ come in fact from the 
already discussed Ref. [20] (Cfr. ^2.l!T ). The analysis of forbidden Xf — > 



final states was the same as for the FCNC ones. Invariant mass distributions of some 



of the analysed states are shown in Fig. 12 where it can be seen that no signal was 
observed. The upper limits on the branching fractions are all of the order of 1CT 5 . 



6 Conclusion 

We have given a broad overview of the current experimental status of searches for 
charm rare decays. While in the past less interest was devoted to these decays than to 
the B and K correspondent ones, a renewed effort is going on both on the theoretical 
and experimental side. This is partially motivated by the recent measurement in the 
other charm processes but also due to new experimental searches going on at LHC 
and B-factories. 

The most stringent limits are still not approaching the Standard Model predic- 
tions and large theoretical errors still dominate the calculations of SM long distance 
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Figure 12: Invariant mass distributions of charged charm hadron decays into SM 
forbidden final states 
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Figure 13: Invariant mass distributions for candidates of D° — > Xii LFV decays from 
the E791 [15] and CLEOII [17J experiments. 

contributions. Nevertheless these searches are already constraining the parameter 
space of various New Physics models. 

Furthermore, as we have seen, there is a large playground of searches, especially 
in the neutral charm sector, which are more than ten years old and could be repeated 
and updated in current experiments, improving the knowledge of this interesting field. 
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